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ABSTRACT 

The ohmic decay of magnetic fields confined within the crust of neutron 
stars is considered by incorporating both the effect of neutron star cooling and 
the effect of space-time curvature produced by the intense gravitational field 
of the star. For this purpose a stationary and static gravitational field has 
been considered with the standard as well as the accelerated cooling models of 
neutron stars. It is shown that general relativistic effect reduces the magnetic 
field decay rate substantially. At the late stage of evolution when the field 
decay is mainly determined by the impurity-electron scattering, the effect of 
space-time curvature suppresses the role of the impurity content significantly 
and reduces the decay rate by more than an order of magnitude. Even with a 
high impurity content the decay rate is too low to be of observational interest 
if the accelerated cooling model along with the effect of space-time curvature 
is taken into account. It is, therefore, pointed out that if a decrease in the 
magnetic field strength by more than two orders of magnitude from its initial 
value is detected by observation then the existence of quark in the core of the 
neutron star would possibly be ruled out. 

Subject headings: magnetic fields — relativity — stars : neutron 
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1. INTRODUCTION 

The magnetic field evolution in neutron stars has been the subject of much discussion 
over the years both in the observational and in the theoretical context. Recent observations 
of pulsars and their statistical analysis seem to imply that isolated pulsars do not undergo 
significant magnetic field decay during their life times (Wakatsuki et al 1992). The possible 
physical processes through which the magnetic field in the crust of isolated neutron stars 
may decay are ambipolar diffusion, Hall drift, Ohmic dissipation, etc. The reviews of Lamb 
(1991), Chanmugam (1992) and Phinney & Kulkarni (1994) provide a reasonable account 
of our present understanding on the decay of magnetic fields in isolated neutron stars. 

Sang & Chanmugam (1987) demonstrated that the Ohmic decay of the dipole magnetic 
fields is not exponential. Considering magnetic field configurations which are initially 
confined to a small part of the crust and which vanish in the stellar core, several authors 
(Chanmugam & Sang 1989; Urpin & Muslimov 1992; Urpin & Van Riper 1993) found 
that the relatively low electrical conductivity of the crustal matter causes the decay times 
too short to be of observational interest if the impurity concentration is high. On the 
other hand, Pethic & Sahrling (1995) showed that if the field does not vanish in the core 
then the shortest possible decay time becomes about two orders of magnitude larger than 
the characteristic time-scale for decay when the field vanishes in the core. Bhattacharya 
& Datta (1996) have studied the ohmic diffusion of magnetic flux expelled from the 
super-conducting core to the crust owing to spindown-induced field evolution and claimed 
that the result obtained by them is in agreement with the ohmic time-scale that appears to 
be required to explain the observed magnetic field strengths of isolated and binary neutron 
stars in the spindown-induced magnetic flux expulsion scenario. 

Recently Sengupta (1997) demonstrated the significant contribution of space-time 
curvature in reducing the decay rate of the crustal magnetic fields in isolated neutron star 
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by assuming a spherically symmetric stationary gravitational field. However, the cooling of 
neutron star has not been incorporated in that work for the sake of simplicity. 

In this paper the Ohmic dissipation of the crustal magnetic field has been investigated 
by incorporating both the effects of neutron star cooling and the space-time curvature 
produced by the intense gravitational field of the star. The paper is organized in the 
following ways: in Section 2, the relevant equations are presented. In Section 3, the adopted 
neutron star model, the initial configurations of the magnetic field and the electrical 
conductivity have been discussed. The results are outlined in Section 4 and finally in 
Section 5 specific conclusions are made. 



Assuming hydrodynamic motions to be negligible and the anisotropy of the electrical 
conductivity of the crustal material is small, the induction equation in flat space-time can 
be written as: 



2. 



MAGNETIC FIELD EVOLUTION IN CURVED SPACE-TIME 





where a is the electrical conductivity. 



The covariant form of Maxwell's equations are given by: 
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The generalized Ohm's law can be written as : 



= (TgTFvxu* 



(4) 



where F^ are the components of the electro-magnetic field tensor, J M are the components of 
the four-current density, u M are the components of the four velocity of the fluid, are the 
components of space-time metric that describes the background geometry and g = det\g^ u \. 
Here and afterwards Latin indices run over spatial co-ordinates only whereas Greek indices 
run over both time and space co-ordinates. 

If a stationary gravitational field is taken into account, then using equations (|2]), ([3]) 
and (|]) (neglecting the displacement current and taking u % = 0) the corresponding induction 
equation in curved space-time can be derived as: 
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For the description of the background geometry I consider the exterior Schwarzschild 
metric which is given by 



ds 2 = (1 - ^)c 2 dt 2 - (1 - ^y'dr 2 - r 2 (d6 2 + sin 2 ^0 2 ), 



(6) 



where m = MG/c 2 , M being the total gravitational mass of the core. The justification for 
adopting the exterior Schwarzschild metric is provided in Sengupta (1997). Since, the crust 
consists of less than a few percent of the total gravitational mass, M can be regarded as the 
total mass of the star. The non-zero components of the orthonormal tetrad of the local 
Lorentz frame for the Schwarzschild geometry is given by 



(*) 



:i - — r l/ \ y {r) = a - ^) 1/2 ' X U = Vr, A 



r 



(0) 



1/r sin 9. 



If Ft a p\ are the components of the electro-magnetic field tensor in a local Lorentz 
frame, then the components of the electro-magnetic field tensor F^s are defined in the 
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curved space-time through the relation: 



F( a 0) - Xj^Xf^F^s- 



(7) 



Using the metric given in equation (0), equation (^) can be reduced to 
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Following the convention, I consider the decay of a dipolar magnetic field which has 
axial symmetry so that the vector potential A may be written as (0, 0, A^) in spherical 
polar co-ordinates where = A(r,9,t). Since the hydrodynamic motion is negligible so 
u l = dx % /ds = and the metric gives 

2m, 



u ' 



1 - 
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Therefore, from equation (|j) we obtain 



dF r(j) c 2 d 



dt 4:7i dr 

and 

dF ff , 



1 , 2m,, /9 „ [ d T. 2m . 1 
-1 1/2 sin 9 — 1 — -F r , 



a 



[ dr 



c 2 d 



dt 4tt 39 



1(1 



a 



2 sinJ# 
or 



r sin^ 



2m x 1 
r sm 9 



+ 







1 



09 \r 2 sin 



9 



1 



(90 \r 2 sin 9 



(9) 



(10) 



Using the definition 

F a /3 = Ap^ a — A a< p 

the above two equations can be reduced to 

'l_!l^l/2 sin 
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for the flat space-time and 

A (j) = -g{r,t) sin 2 6 (13) 
for the curved space-time, where r and 9 are the spherical radius and polar angle respectively 



one gets from equation (P and equation flllj 
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(15) 



respectively, where x = r / R, R is the radius of the star and y = 2m/ R. 



For both relativistic and non-relativistic cases, I impose the usual boundary conditions 
as given in Urpin & Muslimov (1992). 



3. THE MODEL 

I shall follow Urpin & Van Riper (1993) in my approach. The formalisms adopted 
in the present work have been discussed in details by Urpin & Muslimov (1992). The 
gravitational mass of the star and its radius are taken to be 1.4M Q and 7.35 km respectively. 



3.1. The Initial Configurations 

If one assumes the initial value of f(r,t) = f(r) at t = for flat space-time, then for 
curved space-time ([Wasserman fc Shapiro 1983"t Sengupta 1995) 

g (r,0) = g(r) = [r 2 ln(l - — ) +2mr + m 2 }. (16) 

This is because of the fact that any given magnetic field configuration in flat space-time is 
modified by the curvature of space-time produced by the gravitational field of the central 
object. Asymptotically at a large distance g(r) coincides with /(r). 
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I have considered the decay of the magnetic field which initially occupies the surface 
layers of the crust up to a depth (i) x = 0.979 which corresponds to the density 4 x 10 11 gm 
cm~ 3 and (ii) x = 0.9834 which corresponds to density 2 x 10 11 gm cm -3 . 



The electrical conductivity within the crust has been calculated following the 
approaches of Urpin & Van Riper (1993). The net conductivity of the crustal material at a 
given depth is computed as 



where a p h is the conductivity due to electron-phonon scattering and a imp is the conductivity 
due to electron-impurity scattering. The effect of electron-ion scattering has been neglected 
since the region where this effect could be important is sufficiently thin. a imp is inversely 
proportional to the impurity parameter Q defined as 



where n and Z are the number density and electric charge of background ions in the 
crust lattice without impurity Zi and rii are the charge and density of the i-th impurity 
species. The summation is extended over all species of impurity. It is worth mentioning 
that the value of the impurity parameter Q for neutron star crust is not known at 
present. Electron-impurity scattering becomes more important with increasing density 
and decreasing temperature. Hence at the late stage of evolution when the temperature 
becomes low, the conductivity is dominated by electron-impurity scattering. The impurity 
parameter Q has been taken as 0.0, 0.005, 0.05 in the present work. The conductivity is 
calculated by assuming that the region under consideration is isothermal. This is justified 
because the thermal evolution calculations show that the crust becomes isothermal at an 



3.2. The Electrical Conductivity 




(17) 
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(18) 
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age between 1 and a 1000 yr. For the thermal evolution of the neutron star I adopt both 
the standard as well as the accelerated cooling models for neutron stars with gravitational 
mass 1.4M and radius 7.35 km as used by Urpin and Van Riper (1993). The models are 
based on soft equation of state (Baym, Pethick & Sutherland 1971) with normal npe— 
matter in the core and with standard neutrino emissivities for the slowly cooling case and 
with enhanced neutrino emission due to the Urea process on percolating quarks for the 
rapidly cooling case. 



4. RESULTS AND DISCUSSIONS 



Equations fll4|) and fll5 ) with their corresponding boundary conditions have been 



solved numerically using the standard Crank-Nicholson differencing scheme. The results 
for Schwarzschild geometry have been transformed into a local Lorentz frame by using 
equation (0) and are compared with those for the flat space-time in order to visualize the 
effect of general relativity. 

The evolution of the surface magnetic field normalized to its initial value for both 
the general relativistic and the non-relativistic cases with the standard cooling model and 
with the impurity parameter Q = are presented in Figure 1. In this case, the electrical 
conductivity within the crust is determined by electron-phonon scattering only. 

EDITOR: PLACE FIGURE HERE. 

It is shown by Muslimov & Urpin (1992) that the field behavior is qualitatively 
independent of the forms of the initial configurations but the numerical results differ for 
various choices of the initial depth penetrated by the magnetic field. At a very early stage 
of evolution when the crustal matter is melted in the layers of a maximal current density, 
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the magnetic field does not decay appreciably. After the outer crust solidifies, significant 
decay takes place and after 1 Myr no decay occurs if the impurity content is zero. 

The above scenario for flat space-time is not altered with the inclusion of general 
relativistic effects but significant decrease in the numerical value is found when the effect of 
space-time curvature is incorporated. Figure 1 shows that if the initial field configuration 
extends up to the density po = 4 x 10 11 gm cm -3 , then the decay is much slower than that 
for po — 2 x 10 11 gm cm -3 . This is due to the fact that in the former case, the currents were 
localized in deeper layers with higher conductivity from the very initial stage. Figure 1, 2 
and 3 show that during the period when the electrical conductivity is mainly determined 
by the phonon-electron scattering, the decay rate, for the general relativistic case, of the 
magnetic field initially penetrated up to a density p = 2 x 10 11 gm cm -3 is almost the same 
to that when the general relativistic effect is not included and the initial field is penetrated 
up to a higher density po = 4 x 10 11 gm cm -3 . This implies that the curvature of space-time 
acts as if it pushes the currents towards a deeper layer of the crust. 

EDITOR: PLACE FIGURE | HERE. 

EDITOR: PLACE FIGURE | HERE. 

Figure 2 and Figure 3 show the magnetic field decay with the standard cooling 
model and with the impurity parameter Q = 0.005 and Q = 0.05 respectively. The 
impurity-electron scattering is dominant at the late stage of evolution when the crustal 
temperature is low. As a consequence the decay rate after t > 10 Myr changes appreciably. 
Figure 2 and Figure 3 show that the general relativistic effect tends to suppress the role 
of the impurity content that accelerate the decay rate. If the impurity parameter Q is as 
high as 0.05, then at the late stage the reduction in decay rate by the effect of space-time 
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curvature is more than an order of magnitude. Thus, general relativistic effect is more 
important at the late stage of evolution of the crustal magnetic field. 

Figure 4 and Figure 5 present the magnetic field decay with the accelerated cooling 
model and with the impurity parameter Q = 0.005 and Q = 0.05 respectively. If the 
accelerated cooling model is adopted for the description of the thermal evolution of the 
neutron star, then the decay rate even for the flat space-time is found to be negligible 
when the value of Q is zero or very small, i.e., when the electron-phonon scattering is the 
only or the dominating process in determining the electrical conductivity within the crust. 
Therefore, the general relativistic effects in this situation is not interesting. 

EDITOR: PLACE FIGURE | HERE. 



EDITOR: PLACE FIGURE | HERE. 



It is clear from the figures that the magnetic field decay is qualitatively different for 
the standard and the accelerated cooling models. During the neutrino cooling period 
the crustal temperature is much higher for the standard cooling model than that of the 
accelerated cooling case. This makes the electrical conductivity for the standard slowly 
cooling model lower than that for the accelerated cooling model. As a consequence, the 
decay in the magnetic field is less for the accelerated cooling case compared to that for 
the standard cooling model. Unlike the case for standard cooling model, the decay for the 
accelerated cooling model does not exhibit a flat region at 1 - 10 Myr. For the accelerated 
cooling model, noticeable decay occurs due to the electron-impurity scattering only and 
hence there is practically no decay if the impurity content is very small. However, if the 
value of the impurity parameter Q is greater than say, 10~ 3 then after 10 Myr noticeable 
decay takes place. Figure 4 and Figure 5 demonstrate that general relativistic effect reduces 
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the decay rate nearly by 50% even if the value of Q is very high. This decrease in the decay 
rate is almost constant throughout the subsequent period of evolution. But if the impurity 
parameter Q is as high as 0.05, then there is a rapid increase in the decay rate after 10 9 yr 
for the flat space-time. The same type of increase in the decay rate is not found for the 
general relativistic case as the effect of the impurity content is suppressed by the effect of 
space-time curvature. Consequently, if the impurity content is very high, then after 10 9 
yr, the rate of decay reduces by more than one order of magnitude if one incorporates the 
general relativistic effects. However, from the consideration of the age of isolated neutron 
stars, this very late stage difference in the evolution of the magnetic field may not be 
relevant from the observational point of view. Nevertheless, a 50% reduction in the decay 
rate by the effect of general relativity alone is certainly important when accelerated cooling 
model is considered. 

5. CONCLUSIONS 

It is shown that the general relativistic effect plays a crucial role in determining the 
magnetic field decay within the crust of isolated neutron star and the importance of the 
effect is no less than that of the initial depth penetrated by the magnetic field. At the late 
stage of evolution, the effect of space-time curvature produced by the intense gravitational 
field of the star suppresses the role of the impurity content of the crustal material and thus 
reduces the rate of field decay substantially. If the neutron star cools rapidly due to the 
existence of quark matter in the core then the decay rate becomes much lower than that 
when standard slow cooling is considered. A further reduction in the decay rate by the 
effect of space-time curvature results such a small decay in the magnetic field (during the 
whole life span of the star) that irrespective of the uncertainties in the impurity content, 
existence of quark matter that can give rise to accelerated cooling may possibly be ruled 
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out if more than just two orders of magnitude decrease from the initial field strength is 
observationally inferred. 

The author is thankful to the anonymous referee for constructive comments and useful 
suggestions. Thanks are also due to Bhaskar Datta for some discussion and Somnath 
Bharadwaj for a critical reading of the manuscript. 
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Fig. 1. — The evolution of surface magnetic field normalized to its initial value for flat and 
curved space-times with the impurity parameter Q = 0. Solid lines represent the results for 
curved space-time while broken lines represent that for flat space-time. The numbers near 
the curves indicate the value of po/10 11 gm cm -3 where p is the density corresponding to 
the initial depth penetrated by the magnetic field. The results are obtained by using the 
standard cooling model. 

Fig. 2. — Same as figure 1 but with Q = 0.005. 

Fig. 3. — Same as figure 1 but with Q = 0.05. 

Fig. 4. — Same as figure 2 but with the accelerated cooling model. 

Fig. 5. — Same as figure 3 but with the accelerated cooling model. 
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